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The capture and subsequent inspiral of stellar mass black 
holes on eccentric orbits by central massive black holes, is one 
of the more interesting likely sources of gravitational radiation 
detectable by LISA. We estimate the rate of observable events 
and the associated uncertainties. A moderately favourable 
mass function could provide many detectable bursts each year, 
and a detection of at least one burst per year is very likely 
given our current understanding of the populations in cores 
of normal spiral galaxies. 


I. 

The proposed LISA observatory ( ]l] [| and this vol¬ 
ume) has optimal sensitivity to gravitational radiation in 
the 10 -3 10 -2 Hz range. As such it is well matched to 
the expected frequency of gravitational radiation emit¬ 
ted by compact objects in the last stages of spiral-in to 
massive black holes with total masses of Mbh ~ 10 6 M 0 
(see eg. J3J). Previous papers (5|j^] considered the “likely” 
signal from degenerate compact objects in the cusps of 
normal galaxies, captured by central black holes such as 
are inferred to be present in the Milky Way 0 and the 
nearby dwarf elliptical, M32 [^] |l0[ . 

Here we consider the rate for capture of stellar mass 
black holes expected to be present in the cusps of nor¬ 
mal galaxies, and estimate the detectable rate of events 
by LISA, under both very conservative assumptions and 
with more optimistic estimates of the black hole popula¬ 
tion (see also 00 ). 

There are several large, systematic uncertainties in es¬ 
timating the rate of black hole capture. The proportion 
of galaxies harbouring a central black hole is poorly con¬ 
strained, although theoretical prejudices suggest central 
black holes may be ubiquitous 00- The mass func¬ 
tion of any central black holes is highly uncertain; there 
are strong observational biases on detection of central 
black holes, and few galaxies for which there are strong 
observational constraints as yet [jn|. The mass function 
and total number of stellar mass black holes is also highly 
uncertain jl7|. Stellar mass black holes may form with 
masses as low as 2 M 0 and go up to masses of 50-100 M 0 . 
The range of masses of main sequence stars that form a 
black hole remnant and the resultant mass of the rem¬ 
nant is also uncertain, and binarity and metallicity may 
both affect the remnant mass. There are observational 
biases against detecting the higher mass range of stellar 
mass black holes as they may sustain higher accretion 


rates while in binaries, and be detectable in our galaxy 
for correspondingly shorter time. 

We assume, conservatively, that masses of central black 
holes in galaxies are proportional to the mass of the lu¬ 
minous spheroid of the host galaxy |li|, with a space 
density of ~ 3 x 10 -3 Mpc -3 for the range of black hole 
masses of interest; and that stellar mass black holes are 
formed from main sequence stars with zero-age masses 
greater than 15 M 0 , drawn from a Salpeter mass func¬ 
tion extending to 30 M 0 , with a total number fraction of 
~ 10 -4 black holes per formed star. It is possible, if not 
likely, that this underestimates substantially the stellar 
mass black hole fraction. 

In the presence of a central black hole, the stellar pop¬ 
ulation assumes a self-consistent cusp profile, with the 
stellar density increasing like r~ 3 ' 2 ~ p ^ r \ where p(r) de¬ 
pends on the local relaxation time and the original phase 
space density of the stars |ll] |2l|.p|| . For (sub)populations 
of stars for which p(r) > 0, the time scale for two-body 
relaxation, t r , decreases with r (assuming r > r co u, the 
collision radius, at which two-body scattering becomes 
ineffective, for black holes r co u = rs). 

In the cusp of stars around the central black hole, stars 
undergo diffusion in phase space due to three main pro¬ 
cesses: the perturbation of their orbits from the inhomo¬ 
geneities of the potential, due simply to the finite number 
of stars contributing to the mass; large fractional changes 
in energy or angular momentum due to elastic scatter¬ 
ing by close encounters with other stars; and, inelastic 
mergers with other (extended) stars. Due to the pres¬ 
ence of an event horizon at rs, diffusion in phase space 
is not symmetric, rather there is a sink at small radii, 
and stars with angular momentum below some critical 
value are captured by the central black hole (see [p2[|(] 
for discussion). The critical region of phase space for 
which compact stars are captured is well characterised by 
a loss-cone parametrised by a critical angle 0 c (r); here 
we consider the rate for stellar mass black holes in the 
cusp to enter 9 < 9 C through large angle scattering (ie 
80 9 c (r)). The derivation is discussed in |I||j . 

The rate of scattering, R, into the loss-cone from an 
initial orbit of period P, is given by (Jg| 


R(r) 


N*{r)0 2 c (r) / m* \ 2 

P V, Mbh ) 


( 1 . 1 ) 


Where IV* is the number of scatterers interior to r. For 
a given p(r) we can solve for the profile of each stellar 
component of the cusp and hence find the rate for cap¬ 
ture. Typically the total rate is peaked at some radius 
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comparable to or larger than r c , where 9 c {r c ) = 1, unless 
fcoii r c . It is possible to do a more sophisticated anal¬ 
ysis of angular momentum diffusion in the cusps, but the 
rate estimates are dominated by uncertainties in cusp pa¬ 
rameters, composition of the stellar population and the 
black hole mass, not the details of the scattering. 

Stellar mass black holes may be assumed to trace the 
light into the cusp, in proportion to their global number 
fraction. Inside the main-sequence collision radius, the 
black hole fraction increases, as relaxation remains effec¬ 
tive for the compact stars at small radii. For the stellar 
mass black holes in normal nucleated spirals, p(r) ss 0.7 
at r <j 10 6 rg but may flatten to 0.3 inside 10 4 rg. Ne¬ 
glecting loss by capture, we expect high volume densities 
of black holes inside 0.001 pc, even though the total num¬ 
ber of black holes is small. 

The inferred rate of black hole capture from such a 
cusp is 10 -5 y -1 . By assumption there are only a few 
hundred black holes in the cusp, and such a rate is not 
sustainable. Since the relaxation time for the black hole 
population increases with r, we expect the inner cusp 
to be rapidly depleted of > 90% of its low mass black 
holes, with the remainder establishing a steady state pro¬ 
file with capture rate of order 10 -8 y _1 . Typically the 
black hole density profile then breaks at ^ 0.1 pc, with 
the local black hole relaxation time i> 10 9 years. This 
can provide a sustainable dribbling of stellar mass black 
holes from the outer cusp to the inner cusp, with capture 
dominated by scattering from ~ 10 3 " 5 rs. 

LISA strain sensitivity is ~ 10 -23 between 10 -3 and 
10 -2 Hz. The characteristic amplitude from a low mass 
particle, to* of period P is 

h c ~ 4 x lO- 24 ™*/ Mgl (M ag /lO 6 M 0 ) 2/3 (lO 4 s/P) 2/3 

a/Cjr pc 

( 1 . 2 ) 

For ?n* ~ 5 — 10 M 0 , assuming a space density of cuspy 
galaxies with Mbh ~ 1O 6 M 0 of 0.003 Mpc -3 as before, 
we conservatively expect to see about three captures per 
year out to 3 Gpc, if stellar mass black holes are present 
at all in cusps of normal nucleated galaxies. 

The typical initial eccentricity of a stellar mass black 
hole committed to capture by the central black hole is 
ei ~ 0.9995. By construction, we have considered the 
rate for orbits such that dynamical perturbations will not 
cause large fractional changes in the angular momentum 
of the star, after it is scattered into the loss-cone, but 
before it enters the central black hole. The lifetime to 
capture through gravitational radiation emission is well 
approximated by tqw = 10 5 (1 — e 2 ) 7 / 2 (P/10 4 s) 8 / 3 y, or 
~ 10 6 years for a stellar mass black hole with ey = 0.9995 
and apocentre at O(10 3 rg). It is necessary that the pe¬ 
riod shrink to less than 10 4 seconds for the system to 
be detectable by LISA. As the orbit approaches rg, and 
the system departs from the classical regime, the eccen¬ 
tricity is not well defined (see eg. [[□]]). Detailed evolu¬ 
tion of the orbit is complicated, but we can consider the 


qualitative behaviour of the orbit: assuming the peri- 
bothron r p = (1 — e)a for classical semi-major axis a 
and some suitably defined semi-classical eccentricity e. 
Until late stages of the capture, the orbit decay is well 
described by the post-Newtonian quadrupole decay for¬ 
mula, a = — 64/(e)M| ff m*/5a 3 . Then r p = (1 — e)a—ae. 
The first term is negative and much smaller than a, while 
for the orbits we are interested in e < 0 , and thus the 

second term is positive. Hence evolution of peribothron 
will not lead to premature crossing of the event horizon 
by the star. 

In general both the central black hole, and the stel¬ 
lar mass black holes will be rotating and the spins and 
orbital angular momentum of the system will be uncorre¬ 
lated. Given high enough signal-to-noise, either from a 
system with ?n* ~ 5 M 0 and d < 1 Gpc, or to* ~ 50 M 0 
at a redshift of <; 1, post 1 ' 5 and higher order spin-orbit 
couplings may be detected in the frequency change of 
the orbit and the orbit precession with comparable con¬ 
tributions from both the stellar mass black hole spin and 
the central black hole spin (see eg. Q), enabling mea¬ 
surements of general relativistic effects in the strong field 
regime, and correspondingly a unique test of the validity 
of general relativity as a theory of gravity in this regime. 
In addition, any such detection would provide important 
information about the masses and spin of both the cen¬ 
tral massive black hole, and the stellar mass black holes 
present in the centres of galaxies. 

We have assumed that the proportion of stellar mass 
black holes in the high density cores of galaxies is compa¬ 
rable to that estimated for the total stellar population. 
If there is a systematic difference in initial mass func¬ 
tion, theoretical considerations suggest the mass func¬ 
tion should be biased towards more massive stars, and 
hence more massive remnants, in high density environ¬ 
ments [^5|,^6|. An initial higher fraction of black holes 
would lead to more rapid merger (commencing either a 
few million years after the first burst of star formation, if 
the central black hole formed before the stars; or, com¬ 
mencing after the central black hole came to dominate 
the dynamics of the core, if the central black hole arrived 
or formed after the initial stars formed), but the same 
sustainability problems exist independent of the number 
of black holes, provided there are fewer than 10 4 initial 
black holes, and more than 30 or so. Hence most of 
the stellar black holes are swallowed in a short time com¬ 
pared to the lifetime of the galaxy. There are two factors, 
however, which lead to a more favourable rate of bursts: 
with more initial stellar mass black holes, we expect more 
higher (~ 30 — 100 M 0 ) mass ones, and the capture of 
those is detectable by LISA anywhere in the universe; 
it also possible that high mass, low metallicity, stars are 
more likely to form massive black hole remnants. If there 
was even one 50 M 0 black hole in the core of each galaxy 
now containing a ~ 10 6 M 0 central black hole, we ex¬ 
pect several mergers per year at z > 1, all detectable 
by LISA. The redshifting of the gravitational radiation is 
not likely a concern, as central black holes with masses 
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now > 10 6 M 0 must have had somewhat smaller masses 
in the past. 

In our own galaxy, there is now substantial circum¬ 
stantial evidence that there was a recent starburst in the 
central region j27j. If such behaviour is typical of nor¬ 
mal nucleated spirals, then every ~ 10 8 years, the cusp 
population of stellar mass black holes may be partially 
replenished. If this occurs commonly, rates of capture 
in the nearby universe may be an order of magnitude 
greater than we have estimated, and we may expect a 
strong detectable signal every few weeks. 

Other favourable scenarios for providing detectable sig¬ 
nals include the merger of two massive black hole of com¬ 
parable mass, or a massive black hole and an “intermedi¬ 
ate mass” black hole. It is possible that even more mas¬ 
sive black holes form in galaxies ( Mbh ~ 10 1 2 — 10 3 M@), 
either during early stages of galaxy formation, as part of 
the ultra-low metallicity Population III stars, or in clus¬ 
ters of stars, growing from < 10 2 M 0 through mergers 
and accretion. If such black holes can reach the inner 
parsec of their host galaxy, relaxation will bring them to 
small radii and capture by the central black hole. The 
resultant burst of gravitational radiation would be eas¬ 
ily detectable, but both the number of such objects and 
the rate at which they might reach centres of galaxies is 
highly uncertain (see eg. mb¬ 
it seems quite probable that LISA will detect signals 
from the capture of stellar mass compact objects by cen¬ 
tral massive black holes in normal galaxies, if our current 
understanding of the properties of the central stellar pop¬ 
ulation is even marginally correct, and if the masses of 
the central black holes in nearby galaxies are at all typ¬ 
ical. It is possible that the rate of detectable events is 
greater than one per month. The actual detection of 
such events, in addition to testing theories of gravity in 
the strong held regime, would provide the only data on 
the mass function of remnants stars in galactic centres, 
and would probe the low mass end of the range of massive 
central black holes in the distant universe. 
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